NASA Contractor Report 191129 



/ ^ 




An Analytical Optimization of Electric Propulsion 
Orbit Transfer Vehicles 


Steven R. Oleson 
Sverdrup Technology, Inc. 
Lewis Research Center Group 
Brook Park, Ohio 


May 1993 


Prepared for 

Lewis Research Center 

Under Contract NAS3-25266 


NASA 

National Aeronautics and 
Space Administration 


[NASA-CR-191129) AN ANALYTICAL 
DPTIMIZATION OF ELECTRIC PROPULSION 
JRBIT TRANSFER VEHICLES Final 
Report (Sverdrup Technology) 10 P 


N93-27060 


Unci as 


03/20 016A919 




AN ANALYTICAL OPTIMIZATION METHOD FOR 
ELECTRIC PROPULSION ORBIT TRANSFER VEHICLES 

Steven R. Oleson 
Sverdrup Technology, Inc. 

Lewis Research Center Group 
Brook Park, Ohio 44142 

Abstract 

Due to electric propulsion’s inherent propellant mass savings over chemical propulsion, electric 
propulsion orbit transfer vehicles (EPOTVs) are a highly efficient mode of orbit transfer. When 
selecting an electric propulsion device (ion, MPD, or arcjet) and propellant for a particular mission, 
it is preferable to use quick, analytical system optimization methods instead of time intensive 
numerical integration methods. It is also of interest to determine each thruster’s optimal operating 
characteristics for a specific mission. Analytical expressions are derived which determine the 
optimal specific impulse (Isp) for each type of electric thruster to maximize payload fraction for a 
desired thrusting time. These expressions take into account the variation of thruster efficiency with 
specific impulse. Verification of the method is made with representative electric propulsion values 
on a LEO-to-GEO mission. Application of the method to specific missions is discussed. 

INTRODUCTION 

In the coming decades, electric propulsion promises to become a frequently used and highly 
efficient mode of space propulsion. This is due to the inherently higher specific impulse (Isp) of 
electric propulsion (800 to 10,000 s verses 100 to 450 s for chemical propulsion) which produces 
significant propellant mass savings. These savings can then be used to increase deliverable 
payload mass or allow the use of smaller, less expensive expendable launch vehicles. However, 
these savings usually come at the expense of longer mission times. Electric propulsion is already 
being used to provide station keeping propulsion for Earth satellites. In addition, current NASA 
studies have shown the potential for electric propulsion’s use as primary propulsion on 
interplanetary probes, lunar and Mars cargo vehicles, and Mars piloted vehicles. As part of the 
evolutionary process towards the interplanetary application, the initial use of electric propulsion 
(EP) for primary propulsion is likely to be as a highly efficient Earth orbit transfer vehicle (OTV). 

The purpose of this paper is to provide this quick analytical method for predicting EP 
performance, optimal Isp’s, and the corresponding EP vehicle parameters. This method can be 
used to quickly determine the performance of electric propulsion orbit transfer vehicles (EPOTV) 
for any circular orbit to circular orbit transfer mission. By using the predicted optimal Isp method, 
parametric variation of Isp is not necessary thus eliminating one variable from a parametric trade. 

The analysis of EP trajectories is different from chemical impulsive trajectories since the electric 
thrusters operate continuously at very small thrust levels. As a consequence of this, an EP vehicle 
performs a spiral trajectory from orbit to orbit near a planetary body. To accurately analyze this 
spiraling, numerical integration methods are usually used, especially when optimal vehicle steering 
is desired. These methods usually require specific vehicle assumptions which do not allow for 
quick parametric assessments of various thrusters, power systems, power levels, and specific 
impulses. Alternately, by using certain simplifying mission assumptions, analytical methods can 
be used to perform a quick parametric scoping of potential thruster/power combinations, without 
requiring more detailed inputs or specialized user training. In most cases, the efficiency of EP 
thrusters is highly dependent upon the Isp. This variation of efficiency with Isp must be 
considered before choosing thrusters and power levels since it greatly affects payload performance. 
Gilland (1991) did account for this variation for lunar and Mars missions. A similar derivation 
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will be made here for near-planetary EP missions. 

Specifically, this analytical performance technique determines the optimal Isp and maximum 
payload for EP vehicles considering specific thrusters and desired thrusting time. Assumptions are 
presented describing the general characteristics of power source and the thruster systems needed to 
perform a parametric trade. The low thrust mission analysis is simplified to a required low thrust 
velocity increment, AV, to allow a simple analytical relation of the mission requirements to the EP 
vehicle system parameters. EP vehicle parameters are defined and the dependency of optimal Isp, 
optimal power level and trip time on these parameters is then shown. Optimal Isp equations are 
derived for each type of thruster using curve-fit efficiency/Isp relationships and mathematical 
simplifications. The accuracy of these analytically predicted optimal Isp’s is compared to the 
numerically calculated optimal Isp’s using sample value for the mission and vehicle parameters. A 
comparison of the resulting payload fractions is also made. 

ASSUMPTIONS 

The analysis in this paper assumes a single parameter, the specific mass, to characterize the power 
and thruster system scaling. The independent variable is assumed to be the vehicle thrusting time. 
The analysis considers the available power to be constant, which is representative of nuclear 
powered EPOTVs. For solar powered craft, two other affects make the actual trip time longer than 
the thrusting time; (1) the Van Allen belt radiation degrades the power, and (2) shading limits 
power production time to sunlit parts of the orbit. Other methods must be used to assess these 
affects and are not addressed here. 

The efficiency of some electric thrusters is highly dependent upon the effective exhaust velocity, c 
or specific impulse Isp. Depending upon the type of thruster (ion, MPD, or arcjet), the efficiency - 
Isp relation differs. All of the following relations are based upon projected near-term technology. 
The coefficients of these relations are derived from curve fits of empirical data. As long as the 
mathematical form of the efficiency/Isp relationship is valid, however, the method derived in this 
paper is applicable with the appropriate efficiency coefficients. 

For ion thrusters utilizing inert gases (argon, krypton, xenon) the relationship is second-order as 
shown by Equation (1). For MPD thrusters using hydrogen fuel, the relation is the same as (1) for 
Isp’s greater than 7800 s (Gilland, 1991). For Isp’s less than 7800 s the efficiency more closely 
follows relation (2). 

Tlion = (bc2)/( C 2 + d2) [1] T] MPD = ( b C ) / ( C + d ) [2] 

The near-term coefficients b and d for each propellant as well as their feasible Isp ranges are 
presented in Table 1, the corresponding efficiencies predicted by Equation (1) are shown in Figure 
1. Recent H 2 arcjet data (NASA, 1991) indicate that arcjet efficiency does not have a strong 
correlation with Isp. Consequently, it is assumed that H 2 arcjets have a constant efficiency. The 
near-term Isp range and efficiency for H 2 arcjets are also shown in Table 1. 

The nominal range over which an electric thruster’s Isp can be varied is limited. Table 1 shows 
the accepted ranges for each type of thruster. For some orbit transfer missions the optimal Isp may 
be below this Isp range depending upon efficiency, specific mass and desired thrust time. Recent 
research has shown that “derated” xenon and krypton ion thrusters can be run at lower Isp’s with a 
different efficiency relationship than for high Isp operation (Equation 1) (Patterson, 1992). The 
implications of using derated thrusters are not considered in this papa - . 
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MISSION ANALYSIS 


An analytically calculable parameter is used to relate the mission performance requirements to 
vehicle parameters. This mission parameter is the velocity increment required or AV. This AV is 
related to the vehicle performance through the rocket equation: final mass / initial mass = e-^v/c s 
where c is the effective exhaust velocity of the EP thruster. 

Since an EP vehicle thrusts throughout most of its mission, an impulsive AV, (an assumption that 
all thrusting bums change the vehicle’s velocity instantaneously), can not be accurately used. 
However, past work done by Edelbaum (1961) gives an analytical expression which approximates 
the AV required by a low thrust vehicle performing an orbital transfer and a plane change near the 
Earth. This expression is, 

AV = [ V 0 2 + V 2 - 2 V V 0 cos ( jc/2 * e ) ]i/2 (3) 


where V 0 and V are the circular velocities of the original and desired orbits, respectively, and 6 
is the desired plane change. Some of the assumptions accompanying this equation are: (1) the 
orbit must always remain quasi-circular during the transfer, (2) the thrust angle is assumed 
constant during each revolution, and (3) the thrust magnitude is assumed constant throughout the 
transfer. In addition, the thrust- to- weight ratio of the EP system must be 10- 2 or less. 

EP VEHICLE ANALYSTS 

An EP vehicle is modeled by four main elements each of which contributes to the total mass (M c ) 
of the EP vehicle: the propellant, the propellant tankage, the propulsion system, and the payload. 
The following derivation is similar to that done by Stuhlinger (1964). The propellant mass (M p ) is 
the mass of propellant required for the mission. The propellant tankage mass is directly related to 
the mass of the fuel by a tankage factor: K t = mass of tankage / mass of propellant. The 
propulsion system mass (M w ) is defined as the power subsystem mass (main power generation and 
power management and distribution-PMAD) and thruster subsystem mass (including thrusters, 
gimbals, power processing units -PPUs, and propellant feeds but excluding fuel and tankage). 
Finally, the payload mass (Ml) includes useable payload mass, fixed support system masses not 
scalable to the thruster power (for example structures, guidance and control, and communications), 
and any contingency mass. These combine to make the total initial EP vehicle mass: 

hfo ~ h^L + Mp + Kj h4 p> (4) 

Such parameterization of an EP vehicle allows for each component to be compared to the total 
initial mass of the EP vehicle. A system mass fraction (or just system fraction) can be calculated as 
system mass n / total initial mass = JJ. n . The propellant fraction, |i p , is determined easily using the 

rocket equation with the mission’s Edelbaum AV (Equation 3) and the exhaust velocity of the 
electric thrusters (c): 


|i. p = 1 - e- Av / c . 


(5) 


The tankage mass fraction fiKt is: 


P-Kt = Kt Pp . 


( 6 ) 
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The propulsion system fraction |i w is derived by considering the thrust power (Pt): P t = 1/2 ( dm 
/ dt ) c2 , where the mass flow rate is obtained from the propellant required and the given thrusting 
time, x, : dm/dt = Mp / x. The thrust power can be related to the electric power (Pe) input to the 

thrusters by an efficiency factor, q: Pt = T| Pe, where the efficiency accounts only for the actual 
efficiency of the thruster. Line losses, PMAD and PPU efficiencies, and housekeeping power 
requirements are included in the specific mass of the propulsion system. As discussed previously, 
the efficiency of the thruster is known to vary with exhaust velocity or Isp. The specific 
efficiency/Isp relationships for each thruster are described previously. 

Finally, this electric power (Pe) provided to the thrusters is related to the mass of the power 
subsystem by : M w = Pe a, where the term a represents the specific mass of the propulsion 

system as: a = propulsion system mass [kg] / electric power sent to thrusters [kWe]. These 
equations are combined to determine the propulsion system fraction: 

|i w = (a c 2 ) / (2 1| x) (1 - e-Av/c). (7) 

Finally, noting that the values of the system fractions must add to one: M-L + M-w + |i P + P-Kt = 1, the 
payload fraction, fiL, is defined by the other fractions as: 


p L = e-Av/c . [ K t + (a c 2 ) / (2 q x) ] (1 - e-Av/c). (8) 

It is evident that the EP system payload fraction is dependent upon the mission parameters AV and 
specified thrust time x, as well as the system parameters c, a, q, and K t . Consequently, the 
performance of an EP vehicle can be easily predicted for various mission and system assumptions. 

ANALYTICALLY OPTIMAL it L 

Now that an analytical expression for the payload fraction is available, the various relations for 
efficiency (Equations (1) and (2)) are used to find a thruster specific payload fraction: 

JJ-L = e-Av/c . [ K t + (a c 2 ) / (2 q x) ] (1 - e-Av/c) [ q = constant] (9) 

Hydrogen Arcjets 

JJ-L = e-Av/c - [ Kt + a (c 2 + d 2 ) / (2 b x) ] (1 - e-Av/c) [q =( b c 2 ) / (c 2 + d 2 )] (10) 

Ion and MPD Isp>7800s 

Pl = e-Av/c - [ K t + a (c 2 + c d ) / (2 b x) ] (1 - e-Av/c) [q = (bc)/(c + d)] (11) 

MPD Isp< 7800s 

The payload fraction for specific choices of AV, thrust time, specific impulse, tankage, propulsion 
system specific mass, and thruster efficiency factors b and d are calculated using these equations. 
With the choice of a specific mission and a specific system only the thrust time and the constant 
exhaust velocity (or specific impulse) of the thrusters can be varied to alter payload performance. 

Payload optimizations are presented by various authors including Stuhlinger (1964). The op timal 
Isp for maximum payload fraction is shown to be the same as the Isp for the shortest trip time and 
a fixed payload fraction. Stuhlinger developed an approximation of optimal exhaust velocity for 
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ma-rimnm payload fraction as: c opt = ( x / a )V2. However, this parameter does not address the 
dependency of EP thruster efficiency on specific impulse. Gilland (1991) takes these variable 
efficiency equations into account for Mars and Lunar mission scoping. A similar derivation is 
presented here for near-planetary applications. The figures of merit that this method addresses are 
payload fraction and the flight time since both are directly linked to the economics of any mission. 

Equations (9)-(l 1) are easily optimized in terms of exhaust velocity by taking the partial derivative 
of the payload fraction with respect to the exhaust velocity and setting it equal to zero: op.L / o c = 
0. Equation (10) is addressed first since the solution for optimal exhaust velocity or Isp should be 
the same as the optimal for Equation (9) if the efficiency is constant and the efficiency factor d is 
set to zero. The relation giving the optimal exhaust velocity for maximum payload fraction is, 


e Av/c = [Av/c 2 ((bx/a)(l+K t ) + d 2 /2) + c + Av/2] ; (12) 

which can only be solved numerically in its present form Consequently, the substitution, x = Av 
/ c, and the truncated Taylor series: e x = 1 + x + x 2 / 2! + x^ / 3! , are used to simplify Equation 
(12) so that the optimal exhaust velocity, c opt , is the real root of a quadratic: 

c opt = -Av / 2 + [ 2 b x ( 1 +K t ) / a + d 2 -Av 2 /12]i/2. (13) 

Equation (13) gives the optimal exhaust velocity for thrusters having the efficiency relationship 
shown in Equation (1) (current ion thrusters and MPD thrusters with Isp > 7800 s). As a 
consequence of the truncated Taylor series. Equation (13) is not exact. (However, for EP thruster 
Isp ranges of 800 s to 10,000 s and mission AVs below 10 km/s the method’s error is negligible.) 
The approximation loses accuracy at low exhaust velocities and large mission AV’s and thus 
should be checked before using this method. 

By letting thruster efficiency factors d = 0 and b = T) in (13), the optimal exhaust velocity for a 
constant efficiency is expressed as: 

c op t = -Av / 2 + [2 ii x (1 +K t )/a- Av 2 / 12]!/ 2 , (14) 

which may be used for arcjets. Finally, for the efficiency relationship Equation (2), the optimal 
exhaust velocity is: 

c opt = -Av / 2 + [ 2 b x (1+K t )/ct-d Av / 2 + Av 2 / 4 ] t/ 2 , (15) 

which may be used for MPD thrusters with Isp < 7800 s. Due to algebraic complexity, only a 
second order expansion of e x was used, thus making this approximation slightly less accurate. 

VERIFICATION OF THE PREDICTED OPTIMAL ISP 

The accuracy of the predicted optimal exhaust velocity relation for ion propulsion Equation (13) 
as compared to the numerically calculated and the Stuhlinger predicted optimal specific impulses is 
now presented. The accuracy of the other two optimal equations (Equations 14 and 15) is similar. 
The sample mission is a transfer from a low Earth parking orbit to geostationary orbit. It is 
assumed that the initial parking orbit is a circular Earth orbit at an altitude of 500 km with an 
inclination of 28 .7°. The desired geostationary orbit altitude is 35683 km with a 0° inclination. 
For this orbit transfer, the Edelbaum equation produces a required low thrust AV of 5.86 km/s. 
The thruster used to verify the method’s optimal Isp is xenon ion. Xenon ion efficiencies and Isp 
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Toonf e i! n Tabl ® L A mnka § e and reserves fraction of 10% is assumed for xenon 
tw??D 0 i^ r SyStCm SpeClflC I ^ SS is assumed to be 4 0 kg/kWe based upon 
pS P nn^ an \ (Truscello > ct al, 1992) plus an additional 4 
kg/kWe for the PMAD (Gilland, George 1992) A sample thruster specific mass of 6 kg/kWe 

SSn^'r PPU i£^ a 11 ^ P ro P ellant feeds) is added to the power system specific 
masses (Gilland, George 1992). Thus the sample propulsion system specific mass is 50 kg/kWe. 

For the sample mission and system, the optimal Isp’s predicted by Equation (13) are shown in 
Figure 2 as a function of thrusting time. The optimum Isps predicted by Stuhlinger and the 
numerically calculated optimum Isps are also shown. The equations derived in this paper show 

° PtimUm SpeciflC in * uIse - Stuh ^er’s Isp varies 

IJe corresponding maximum payload fractions found using the predicted, Stuhlinger, and actual 
optimum Isp s are shown in Figure 3. It is clear, at least for the sample case chosen, that the 
maximum P^oad fraction predicted by the analytical method is very close to the actual optimum 
payload fraction. While the maximum payload fraction predicted by Stuhlinger shows only a small 
misleadhfg ^ rcSultUlg s y stem parameters, such as required power and amount of propellant, are 

APPLICATION OF THE METHOD 

The payload fractions for the sample mission using the three thruster types are shown together in 
Figure 4. Figure 5 shows corresponding optimal Isp’s. By using these figures, a preliminary 
scoping of the mission performance of each thruster type is available. For EP systems a trade 
exists between mission duration and payload delivered. While no ‘optimal’ payload and nip time 
is evident, designing a vehicle with attributes from the ‘knee' of the curve, (where the ret^rf 
more payload for longer trip time diminishes) is suggested. In further analyzing the mission a cost 
analysis should be done balancing the operations cost against reduced launch cits. 

It is evident that for short thrusting times, some of the optimal Isp’s are below the thruster Isd 
range provided with Table (1). This indicates that either the lowest Isp boundary value must be 

ISE— must be done t0 pemut °P eration in this Is P range, if the performance advantage 
CONCLUSIONS 

A quick analytical method for predicting mission performance for EPOTV vehicles is shown Rv 
using the optimal exhaust velocity (Isp) equations derived in this paper, one p^ter of defining 
the mission (Isp) need not be parametrically varied. The optimal Isp equations are shown to 

accurately 6 ° ptimal Isp for maximum payload fraction for a desired thrust time quite 
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Nomenclature 


English 

c: Exhaust Velocity 

go: Reference Gravitational Acceleration 

Isp: Specific Impulse, c / go 

b: Thruster variable efficiency factor 

d: Thruster variable efficiency factor 

V: Circular Orbit Velocity 

K t : Tankage Fraction 

M: Mass 

P t : Thrust Power, 1/2 mass flow rate * c 2 

dm/dt: Propellant Mass Flow Rate 
P e : Electric Power Produced by Power 

Subsystem 


Greek 

T) : Propulsion system efficiency, P t / P e 

a: Propulsion (Thruster and Power) 

System Specific Mass, M w / P e 

0: Inclination change (*) 

x : Thrusting Time 

(i: Mass Fractions, subsystem Mass / Mq 

AV: Velocity increment 

Subscript 

o: initial 

L: Payload, Support, and Contingency 

w: Propulsion Subsystem 

p: Propellant 

opt: Optimal 


7 


Payload Fraction Thruster Efficiency, 


j r.i i rm n iTjrrrrra « i >fqnvwn»i.f. ■ f-™ ?*ttt 1 hhhiim— « 

Thruster 

Type 

Efficiency 
Factor b 

Efficiency 
Factor d (m/s) 

Nominal Isp 
Operation Range 

"Derated" Isp 
Operation Range 

Xenon Ion 

0.81 

13500 

3500 - 5000 s 

800 . 3300 s 

Krypton Ion 

0.81 

17000 

5000 - 8000 s 

1500 - 5100 s 

Argon Ion 

0.81 

24500 

7000 - 10000 s 

5000 - 7000 s 

H2 MPD: Isp < 7S00 s 

0.86 

25900 

2000 - 7800 s 


H2 MPD: Isp > 7800 s 

0.92 

50300 

7800 - 10000 s 


Hydrogen Arcjet 

Constant efficiency * .33 

1000 - 1500 s 




Isp (s) 

FIGURE X. Various Thruster Efficiencies 

Predicted by Curve-Fit Relationships. 



FIGURE 2. Comparison of Predicted, Stuhlinger 

and Actual Optimal Isps for a Xenon Ion Thruster 



FIGURE 3. Comparison of Predicted, Stuhlinger, 

and Actual Optimal Isps on Payload Fraction 
for a Xenon Ion Thruster 



Thrusting Time (days) 

FIGURE 4. Maximum Payload Fraction verses 

Thrusting Time For a LEO-to-GEO Mission 



Thrusting Time (days) 

FIGURE 5. Optimal Isp’s vs Thrusting Time 
For a LEO-to-GEO Mission 
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